Numerical solution of a fuzzy time-optimal control problem 
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Abstract. In this paper, we consider a time-optimal control problem 
with uncertainties. Dynamics of controlled object is expressed by crisp linear 
system of differential equations with fuzzy initial and final states. We introduce 
a notion of fuzzy optimal time and reduce its calculation to two crisp optimal 
control problems. We examine the proposed approach on an example. 

Keywords: Optimal-time control, fuzzy set, maximum principle, mathemati- 
cal pendulum. 

1. Introduction 

Many researchers investigate optimal control problems with uncertainties. In [1], 
Gabasov et al. consider optimal preposterous observation and optimal control prob- 
lems for dynamic systems under uncertainty with use of a priori and current infor- 
mation about the controlled object behavior and uncertainty. In |2j, Gabasov et al. 
investigate for an optimal control problem under uncertainty the positional solutions, 
which are based on the results of inexact measurements of input and output signals 
of controlled object. In jSj, Gabasov et al. study a problem of optimal control of a 
linear dynamical system under set-membership uncertainty. 

Fuzzy time-optimal control problem is investigated in different forms in [l]-[6]. 
In [1], Plotnikov proves necessary maximin and maximax conditions for a control 
problem, when behavior of the object is described by a controllable differential inclu- 
sion with multivalued performance criterion. In f5], Sakawa et al. propose a fuzzy 
satisficing method for multiobjective linear optimal control problems. To solve these 
problems, they discretize the time and replace the system of differential equations by 
system of difference equations. In [6J , Molchanyuk and Plotnikov study the problem 
of high-speed operation for linear control systems with fuzzy right-hand sides. For 
this problem, they introduce the notion of optimal solution and establish necessary 
and sufficient conditions of optimality in the form of the maximum principle. 
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In this paper, we consider a time-optimal control problem with crisp dynamics 
and with fuzzy start and target states. We interpret the optimal time as a fuzzy 
variable and propose a numerical method to calculate it. 

The paper consists of 5 sections. In Section 2, we describe the classical time- 
optimal control problem. In Section 3, we define the fuzzy time-optimal control 
problem and propose a method for calculation of fuzzy optimal time. In Section 4, 
we show the proposed approach by an example. Finally, we give concluding remarks 
in Section 5. 

2. Classical linear time-optimal control problem 
Let the behavior of a controlled object is definite (crisp) and described by the following 
linear system of differential equations: 



Here x is n-dimensional vector-function that describes the phase state of the object, 
A is an n X n matrix, u is n-dimensional control vector-function. 

Let U C i?" be a nonempty compact set. If measurable function m, defined on the 
interval / = [toi^i], satisfies the condition u{t) G U for each t & I, then u is called as 
admissible control. It is known that for any admissible function u and for any initial 
state p the initial value problem 



has a unique solution [?]. This solution x describes how the phase state changes under 
the influence of admissible control u. 

Assume that the start time to and the start state p are given. If we want to 
transfer the object to a given state q in the shortest time by choosing an appropriate 
admissible control u, we have the following Classical time-optimal control problem of 
1st type: 



X = Ax + u 



(1) 



X = Ax + u 



x{to) = p 



ti — to ^ min 



(2) 



u 



Subject to 



X = Ax + u 





(3) 
(4) 
(5) 



Note, that the flnish time ti is not known beforehand and is determined as a result of 
solving the problem. Summarizing, 1st type Classical time-optimal problem (I2])-(|5]) 
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is a problem of finding an admissible control u, which transfers the system from the 
initial phase state p to the final phase state q in the shortest time. 

Now, let nonempty compact sets Mq and Mi from i?", an interval / = [to;^i] ? and 
an admissible function u on this interval are given. If the system ([T]) has a solution 
x{t) such that x(to) £ Md and x{ti) e Mi, then it is said that the control function 
u transfers the object from the initial phase set Mq to the final phase set Mi on the 
interval [to, ^i]- If we want to transfer the object from the set Mq to the set Mi in the 
shortest time, we have the following Classical time-optimal control problem of 2nd 
type: 

^1 — ^0 — ^ niin (6) 

u 

Subject to 

X = Ax + u (7) 

x{to) e Mo (8) 

xih) G Ml (9) 

where Mo and Mi are given start and target sets. The solution u of the problem 
([n])-® is called optimal control. The solution x of the system ([7])- ([21), corresponding 
to the optimal control u, is called optimal trajectory. If u{t) is an optimal control and 
x{t) is a corresponding optimal trajectory, then {u{t),x(t)) is called to be an optimal 
pair. 

We note that the Classical problem of 2nd type can also be reformulated as follows: 



ti - to ^ min (10) 

X = Ax + u (11) 

x{to)=p (12) 

x{h)=q (13) 



2nd type Classical time-optimal problem ([H])-© (or f[T(l - f[T^ ) is well studied [7]. 
Below we give necessary conditions of optimality for this problem [Tj. 

Definition 1. (Maximum principle). Let u he an admissible control defined on an 
interval [to,ti] and let x be a solution of the system (13)- We say that the pair 
{u{t) , x{t)) satisfies maximum principle on the interval [to, ti] if the conjugate system 

has such a nontrivial solution ip = {ipi, ^"2, V'n) that the following conditions hold: 

1) maximum condition: {u{t)^ip{t)) = c{U,ip{t)) for almost any t G [to,^i]; 

2) transversality condition on Mq: {x(tQ),ip(to)) = c(Mo, ^/'(to)); 

3) transversality condition on Mi: {x(ti), —ip(ti)) = c(Mi, —ipiti)). 
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Here (m, ■?/') = uiipi + U21P2 + ••• + Unipn denotes the inner product of vectors u and 
ip from i?" and c{S, ip) = max (s, ip) denotes the support function of the compact set 

S from i?". 

Theorem 1. ^^7] (Necessary conditions of optiniality for the time-optimal control 
problem). Let Mq and Mi be nonempty convex compact sets. Also let the function 
u defined on [to,^i] be an optimal control for the problem and x be a corre- 

sponding optimal trajectory. Then the pair {u{t),x{t)) satisfies maximum principle 
on the interval [to;^i]- 

3. Fuzzy linear time-optimal control problem 
If the start and target values in Classical problem of 1st or 2nd type are fuzzy, we 
obtain the following Fuzzy time-optimal control problem: 

^1 - ^0 min (14) 

u 

Subject to 

X = Ax + u (15) 
x{t,)=l (16) 

x(ti)=C (17) 

where ^ and C, are given fuzzy initial and final vectors (or sets). 

Depending on different definitions for derivative of fuzzy function or different 
definitions for solution of system of differential equations, the problem flT4|) -f lT7|l can 
be interpreted by different ways. We will interpret the problem f lT^ -f lT7|) as a set of 
1st type Classical problems ([2])- ([5]). Each problem is obtained by taking the initial 
value p from ^ and the final value q from C. We denote by ti pq, Upg and Xpg the 
solutions of the problem ([2])-([S]). Let a = min |/i^(p), /i^(g)} (where n^{p) denotes 
the membership of p in ^). We call itipg,Upg,Xpg) to be a solution of the problem 
([T1D-([I7D with possibility a. 

Set of all ti^pg, defined above, determines a fuzzy number ti. We will investigate 
how to calculate ti. Functions ti(a) and ti(a), which indicate the left and right 
boundaries of a-cuts, determine the number ti fully. Thus, the problem of calculation 
of fuzzy optimal time is reduced to calculation of the functions ti{a) and ti{a). 

As it is known, the initial and final values of the optimal solution x{t) of the 
problem (|n])-(in]) are achieved on boundaries of the sets Mq and Mi [7]. Thus, value 
of ti{oi) can be obtained by solving the problem ([6])- ([9]) with taking Mq = ^„ and 
Ml = (^^ and denote a-cuts of ^ and (, respectively), namely the problem: 



^1 — ^0 — ^ min 

u 



(18) 
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x = Ax + u (19) 

^(^o) e e„ (20) 
e C (21) 

The problem ( |T8l) -( l2Tl) is a classical problem of 2nd type. 

Taking into account f ITU]) . it can be seen that ti(a) = min ^i,pg- Note that, 

the value ti (a) means the shortest time between two points, one of them is from the 
set and another is from in the best case. Similarly, ti{a) means the shortest 
time in the worst case: 

ti{a) = max ti^pg (22) 

Calculation of ti{a) can be performed similarly to calculation of ti{a). 

4. Example 

In this section, we apply the proposed approach to a fuzzy time-optimal control 
problem. The problem is a fuzzified version of the crisp problem of damping of 
mathematical pendulum, presented in [7j. 

Example 1. Solve the fuzzy time-optimal control problem (Note that below = 0): 
ti — )■ min 

u 

Xi = X2 

X2 = -Xi + U2 

U = {U= {ui,U2)\Ui = 0, \U2\ < 1} C i?2 

2^1(0) = ^1 = (-6, -5, -4); X2(0) = ^2 = (2, 3, 4) 
Mh) = Ci = (-0.5, 0, 0.5); X2(ti) = C2 = (-0.5, 0, 0.5) 

Here ^,1,^2X1 ^i-nd C2 S'-i's triangular fuzzy numbers. 

Prom above = {(xi,X2)|a; — 6 < a;i < —4 — a, a + 2 < X2 < 4 — a} C i?^ and 
Ca = {(2^15 3;2)|tt — i < xi < 1 — a, a — 1<X2<1 — «}C R'^. 

System's matrix is A = ^ ^ ^ . Since, —A* = A (here A* denotes the conju- 
gate matrix of A), the conjugate system is: 
^1 = ^2 

Support function of U is c{U,tlj) = \'ip2\- Then, from the maximum condition 
{u{t),ip(t)) = c{U,ip{t)) we have U2{t)ip2it) = |V'2(^)I- Consequently, 
U2{t) = 1, if 7/'2(t) > 0; 
U2{t) = -1, if V'2(i) <0; 
-1 < U2{t) < 1, if 4j2{t) = 0. 
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Figure 1: An optimal trajectory can be realized by combining of clockwise motions 
on circles with centers K and L. 

Let us find the solution of the conjugate system corresponding to an initial condi- 
tion V^(0) G C, where C is the unit circle. Initial point can be represented in the form 
of = (cos a, sin a) with a e [0, 2n). Then the solution of the conjugate system is 
il^iit) = cos{a — t),'ip2{'t) = sin(a — t). The function ip2i't) = sin(Q; — t) changes its sign 
for first time at r < tt (r = a, if < a < tt and r = a; — vr, if 7r<a< 27r) and then 
after each tt time period. Depending on a, the sign of the function il'2{t) = sin(a — t) 
in the interval [0, r] is either positive or negative. Thus, according to the maximum 
condition, the initial value of the optimal control U2{t) is either 1 or —1. After r < vr 
time units, it switches from 1 to —1 or vice versa. Then, it repeatedly changes its 
sign after each tt time period. 

Below we interpret the behavior of the object as a motion of the object in the 
phase plane XiX2- 

Solutions of dynamic system corresponding to U2{t) = 1 are in the form x{t) = 
(1 -|- ccos{ip — t),csin((y9 — t)). In the phase plane these solutions constitute 
concentric circles with center at L(1,0) (Fig. 1). The motion on these circles is 
clockwise with constant speed and whole turn takes 2tt time units. 

Similarly, solutions of dynamic system corresponding to U2{t) = — 1 are in the 
form x{t) = (—1 + ccos{ip — t), csin(<^ — t)). In these solutions constitute circles 
with center at K{—1,0) (Fig. 1). The motion on these circles is clockwise with 
constant speed and whole turn takes 2n time units. 

Note that angular speed is w = 1 for both motions mentioned above. So, the 
angle formed by the object during its motion and passed time are equal in value. 

Let us emphasize two facts which will be used in arguments below. 1) In circular 
motion with a; = 1 after tt time period the object will be in the position which is 
central symmetric point of the previous position. 2) The symmetric point of (a, b) is 
(—a — 2, —b) with respect to center point K. If center is L, then the symmetry of 
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Figure 2: A sample of optimal trajectory with 3 switches. 

point (c, d) is (— c + 2, —d). 

Now we investigate how is a motion of the object corresponding to an optimal 
control in the phase plane for a start point S and a target point T. Let us consider 
the case when the object starts with control u = —1 (The case with start control 
u = 1 can be investigated similarly). Let k denote the number of control switches. 
We consider the cases k = (motion without switch) and k >1 separately. 

In the case k = 0, running from the start position 5" and moving along a circle 
with center K the object reaches the target position T. This case occurs, only if 
\KS\ = \KT\ (Here \AB\ denotes the length of the segment AB). The motion time 
is ti = 6^ = LSKT (Here LSKT denotes the value of the angle SKT). 

Now let A; > 1. We differ the cases when k is odd and when k is even. 

Let us consider the case that k is odd number and take = 3 for clarity. The 
object runs from the point S along a circle with center K and after r time period 
arrives a point Xi{x^y) (Fig. 2). The points S and Xi are on the same circle. 
Consequently: 

\KXi\ = \KS\ (23) 

At the point Xi the control switches for the first time and becomes u = 1. Under 
this control, the object moves along a circle with center L. After tt time units it 
arrives a point X2{—x + 2, —y). Here the control switches for the second time and 
under new control u = —1 (moving on circle with center K) after vr time the object 
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Table 1: Point of k-th control switch for optimal motion 
k (odd) Xk k (even) Xk 
i i^) 2 (-x + 2,-y) 

3 (x-4,?/) 4 (-X + 6, -y) 
5 lx--8,y) 6 {-x + 10,-y) 
7 ix-12,y) 8 i-x + U,-y) 

reaches a point X^. = X-j{x — A,y). At the point X^ the control switches for last time 
and becomes u = 1. The object continues its motion on a circle with center L up to 
the target point T. For the aforementioned motion, the points Xk and T must be on 
the same circle with center L, i.e., 

\LXk\ = \LT\ (24) 

It can be seen from Table 1 that for an odd k (including k = 1) the last point of 
control switch is 

Xk = ixk,yk) = ix-2{k-l),y) (25) 

Let S = ipxiVy) and T = {qx,qy)- To calculate unknown coordinates x and y 
we use equations (12^ and (^^. Using (1231) . these equations can be rewritten in 
coordinates as follows: 

(a; + l)2 + y2 = rl = {px + iy+pl (26) 
{x+l-2kf + y' = rl = {qx + lf + ql (27) 

Subtracting (l27|l from (126|) we have: 4A;(x+ 1) — 4fc^ = rf—r"^. Then we can determine 
X and ?/ as follows: 

X = \, ^ +k-l (28) 
Ak 



y = ±^rl-{x + lY (29) 
If X and ?/ have been determined we can calculate the passed time: 

ti = iSKXi + {k- 1)71 + iXkLT (30) 
Let us find an evaluation for k. From (128|1 and ( 1291) we have 
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4 

Hence, we obtain the following evaluation 

fcmin = r|ri - ral /2l < A; < [(n + r^)l1\ = k (31) 

where \x] and [x\ denote ceiling and floor of x, respectively. By taking k = /Cmin, we 
have a feasible motion. Hence, using formula (l30l) . we get: 

ti,opt <ti = iSKXi + (fcinin - l)vr + LXkLT < 27r + (fcmin - l)vr + 27r 

Then, we have /Cmax < ^min + 4 <^==^ fcmax < ^min + 3. Consequently, we obtain the 
following upper evaluation for A;, by using (13T|) : 

A^max = min "l/^min + 3, k 

The case when A; > 1 and k is even can be investigated by similar way. In this 
case the last point of the control switch is (Table 1): 

Xk = {xk,yk) = {~x + 2{k-l),-y) (32) 

The last control is m = — 1 and, consequently, the object finishes its motion on a circle 
with center K. Hence, r| = {Qx — + Qy- Except this value, the formulas for x and 
y become the same as ( l28l) and ( l29l) . The motion time is: 

ti = iSKXi + {k- l)n + IXkKT (33) 

Above we have investigated the case when the start control u equals to —1. In 
the case where m is 1 we have the following final formulas: 

rl = (p,-l)2+p2 (34) 



(g^. + 1)2 + ql if k is odd 
{qx — 1)2 + qy, if k is even 



(35) 



Y - {x + 2{k-l),y), if Zeis odd , . 

Xu - [x,,yu)-\ (_^_2(A;-l),-y), if A; is even 



X = -i^J^ + k-l] (37) 



y = ±^/rl-ix-iy (38) 

iXkKT, if k is odd 
LXkLT, if k is even 



+ /QTv ^ n 1^ , J ^XkKT, if A; is odd . . 

h = iSLXi + (A; - l)7r + <^ / ;f (39) 
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Figure 3: The membership function of fuzzy optimal time ti. 

The above formulas, given for different situations, were obtained on the base of the 
necessary conditions for optimality. Therefore, every solution constructed on these 
formulas may not be optimal. However, the optimal solution is among all solutions, 
constructed for different start controls and for different values of k. 

Based on the above arguments and formulas a computer program is implemented 
to calculate the optimal control for a given pair of start point 5* and target point 
T. Firstly, by taking start control u = —1, after taking u = 1 and in both cases by 
changing the value of k from kjam to /Cmax a solution is constructed (if there is any). 
The solution with the shortest time is the optimal solution, transferring the object 
from S to T. 

Now, let us describe how we calculate the fuzzy optimal time ti numerically. To 
calculate the value we place equally spaced nodes on the boundaries of the 

regions C,^ and The shortest time among all possible start-destination node pairs 
(p, q) gives the approximate value of ^1(0;). 

To calculate the function we discretize the problem (l22ll and solve it numer- 
ically. 

The membership function of fuzzy optimal time ti, obtained from calculations, 
is depicted in Fig. 3. The value ti ~ 8.78 with possibility 1 corresponds to the 
solution of the crisp problem {p = (—5,3) and q = (0,0)). The least value ti ~ 5.97 
with possibility occurs when p = (—4,2) and q = (—0.5,0.5). The largest value 
ti ^ 11.76 with possibility corresponds to the pair p = (—6,4) and q = (0.5,0.5). 
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5. Conclusion 

In this paper, we investigate the problem of time-optimal control with fuzzy initial 
and final states. We interpret the problem as a set of crisp problems. We perform 
the calculation of fuzzy optimal time by solving crisp optimal control problems of two 
types. We exhibit the proposed approach on a numerical example. 
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